Global Dynamics of Cosmological Expansion with a Minimally Coupled Scalar Field 



> 



David I. Santiago* 
Department of Physics, Stanford University 

Alexander S. Silbergleiv 

Gravity Probe-B, W.W.Hansen Experimental Physics Laboratory, Stanford University 

(February 7, 2008) 

We give a complete description of the asymptotic behavior of a Friedmann-Robertson- Walker 
Universe with "normal" matter and a minimally coupled scalar field. We classify the conditions 
under which the Universe is or is not accelerating. In particular, we show that only two types 
of large time behavior exist: an exponential regime, and a subexponential expansion with the 
logarithmic derivative of the scale factor tending to zero. In the case of the subexponetial expansion 
the Universe accelerates when the scalar field energy density is dominant and the potential behaves 
in a specified manner, or if matter violates the strong energy condition p + 3p > 0. When the 
expansion is exponential the Universe accelerates, and the scalar field energy density is dominant. 
We find that for the Big Bang to occur at zero scale factor, the equation of state of matter needs 
to satisfy certain restrictions at large densities. Similarly, a never ending expansion of the Universe 
constrains the equation of state at small matter densities. 
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In this letter we give a complete generic description of the late time asymptotic behavior of the cosmological 
expansion under quite general conditions. While there is literature devoted to specific potentials and interesting 
attractor solutions within those models in the present work we study the generic features of the late time 

asymptotics for any potential. We choose units in which G — c = 1. 

Recent results from supernovae type la jl| indicate that the Universe may be accelerating. If this is the case, we are 
currently entering a period of cosmological inflation. An accelerating Universe is sometimes interpreted as evidence 
for a cosmological constant, which is not necessarily so. Rather, this testifies only that the dominant material in the 
Universe is characterized by an equation of state that satisfies p + 3p < p . "Exotic" matter described by this type 
of equation of state has been called quintessence |l|] . A self-interacting scalar field is one way to provide acceleration 
without assuming this peculiar property of "normal" matter Jl|-§|. As far back as ten years ago, Ratra and Peebles 
considered a cosmology where the scalar field energy density becomes dominant at the present cosmological epoch 
They found asymptotically stable equilibrium solutions in which the scalar field energy density dominates the 
| | dynamics. Their work has recently been resurrected by Liddle and Scherrer 0. These results correspond to certain 
bJ[)' special forms of the scalar field potential and also rest upon other assumptions. 

Since the Universe is very homogeneous and isotropic on large scales, it can be approximated by a Friedmann- 
• rH , Robertson- Walker metric: 

b. ds 2 = -dt 2 + R 2 {t)[-^—:+r 2 dtt 2 ) , (0.1) 



where dft is the element of solid angle and k = -1, 0, or 1 according to whether the Universe is open, flat, or 
closed. In this work we consider mainly open and flat cosmologies. The only energy-momentum tensor consistent 
with homogeneity and isotropy is the one corresponding to a perfect fluid. Therefore we consider a Universe filled 
with material of energy density p and pressure p, and with a self-interacting scalar field tp that has potential V(ip). 
We assume that 

p > (0.2) 
p + p>0 (0.3) 
V(ip) > (0.4) 



Those are fairly reasonable restrictions. Condition ( p.2| ) states that the energy density is positive, while ( |0.4| ) ensures 



that the scalar field energy density never becomes negative. We believe it is unreasonable to consider material with 



* david@spacetime .Stanford .edu 
t gleit @relgy ro .stanford.edu 



1 



pressure more negative than the vacuum (p = —p) and therefore we impose ( |0.3| ). We also exclude the vacuum case 
itself because the cosmological constant can be treated as part of the scalar field potential. The matter p in the 
Universe is made up of baryons, photons, neutrinos, dark matter, and whatever else it might be. We suppose the 
pressure to be some function of the density, p = p(p). The equation of state in the present epoch is very close to that 



of a pressureless fluid ("dust"), but we consider an entirely general equation of state, assuming only (0.3) and 



p(0) = (0.5) 



(in the absence of matter there should be no pressure at all). 
The field equations are 

^=8irp + ip 2 + V(ip)-3^^p T -3^ (0.6) 

P=-^§(P+P) (0.8) 

where the total energy density px — 8np + p v , and p v = ip 2 + V((p) is the energy density of the scalar field. The 
equations also combine to give 



-4tt (p + 3p) - - (p v + 3p v ) , (0.9) 



where p v = ip 2 — V(ip) is the scalar field pressure. Thus the expansion is accelerating when (p v + 2>p v ) < and larger 
in magnitude than 87r (p + 3p) . 

We consider an expanding (as opposed to contracting) Universe. Hence we choose the positive square root in 



equation (|0.6|), which, in view of equations (|0.2j), ( |0.4| ) an d k = 0,-1, implies R/R > 0; hence R > 0, so R — > oo as 
t — ► oo. We thus have a 4D dynamical system whose trajectories {R(t), p(t), (p(t), ip(t)} are specified by initial values 
R > 0, p > 0, ipo, ip at some t = t . 

The energy conservation equation fl0.8| ) can be integrated to 

d ^ 3i n r#) (o.io) 



By equations ( JO. 8|) and ( |0.3| ) , p(t) is a monotonically decreasing (positive) function, therefore it has a nonegative limit 
when t — > oo. Ho weve r, for l arge times R — > oo, so the integral in the left-hand side of ( p. 10 ) must diverge. Because 



of our conditions (3.3) and ( p.5|) , this can only happen if limt^oo p(t) = 0; the matter density in the Universe goes 
monotonically to zero with its expansion. Note that the divergence of the integral (0.10) at p — requires the pressure 
go to zero fast enough for small densities (c. f. the integral converges if p oc p 7 , 7 < 1 as p — * +0). It is sufficient for 
the divergence that the derivative p'(0) exists, which we assume in the sequel. 

We now briefly turn to the Big Bang. We require that R — > at the time (either finite or —00) of the BB singularity 



in the past. As Rq — > 0, In (R/Rq) — > 00, and again the integral in ( 3.10 ) must diverge, which is impossible unless 
p — > 00 (a definition of the BB singularity). However, this integral does not diverge at po — 00 for any thinkable 
dependence p(p): for the divergence, the pressure must not grow too fast at large densities. Thus a surprising fact is 
that the very existence of the BB imposes a limitation on the equation of state in the early Universe. In terms of a 
power scale, if p oc p 1 for p ^> 1, then 7 < 1. It is suggestive that if 7 > 1 not only does the integral converge, but 
causality is violated: the speed of sound is 

2 f dp\ x 



dp 

so it becomes greater than one for sufficiently large p. 

Once again, for the BB to occur at Rq — and/or the infinite expansion to exist, the equation of state must ensure 



that the integral (0.10) diverges at po = 00 and/or p = 0. This result does not depend on the scalar field system, as 



it is minimally coupled; in particular, it is valid for V = as well. 
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While certainly thought provoking, these considerations about the BB should be taken with caution as the situation 
in the very early universe is probably considerably more complicated. In fact, at very early times the universe probably 
underwent a period of inflation |^ [tJ. In such a period the evolution of the universe is not dominated by the energy 
density of regular matter, but by the energy density of the inflat ion field(s). In fact inflationary scenarios in which 
the expansion is exponential violate the weak energy condition (0.3). Since inflation is the only known solution to 
the horizon problem our assumptions are probably too simplistic when applied to an early universe which 

underwent an exponential inflationary phase. On the other hand, acce lerated sub-exponential expansion that solves 
the horizon problem can occur in the early universe and condition (3.3) can be satisfied. Global constraints regarding 
the weak energy condition and generic inflationary scenarios have recently been studied by Vachaspati and Trodden 
H .Whether the existence of the BB at Rq = implies an equation of state which is causal is a matter of a separate 
investigation. 

In order to study the limit of large times for the cosmological dynamical system (0.6) — (3.8), we first note that 
equations (|0.8|) and (p.7[) combine into 



PT = [4tt (p + p) + tp 2 } < 0, 



(0.11) 



with the inequality implied by ( |0.3| ) (it is straightforward to see that the density does not turn to zero at any finite 
time). This shows that the initial store of energy is dissipated into the expansion, and that, just as p, the total 
energy density px is a (positive) decreasing function of time; it is a Lyapunov function of our dynamical system (e. 
g. ref ||, 2.3). Thus pr has a nonegative limit lim^oo pt — p^ = p^(Ro, po, <Po> <Po) > 0; since p — * 0, we have also 
lim^oo p v = p%? . Therefore the large time behavior of the cosmological expansion is described by 



E>2 

p^0 

Pip = p 2 + V(y) =/>§? + ..., 
where dots stand for the terms tending to zero; in addition, equations ( |0.7D and (0.9) become 

1 dV{tp) 



P + f3p + 



2 dip 



+ ... = 0, f3=^pf>Q 



(0.12) 

(0.13) 
(0.14) 



(0.15) 



3- = Pt ° - 3tp 2 + . . . = W(ip) - 2pf? + . . . (0.16) 
R 

Depending on whether pj? is positive, there are clearly two different large time regimes of the expansion. 

If Pt > 0, the Universe expands (and of course accelerates) exponentially; this is the only possible regime when 
the potential is bounded away from zero, V(p) > V m i n > 0, for instance, if a cosmological constant exists. 

The expansion is subexponential with R/R — > when p™ = 0, which requires the zero value to belong to the 
closure of the V((p) image. In this case, the Universe could accelerate or decelerate depending on further properties 
of the scalar field potential. 

We first consider the subexponential case, i. e., any solution with p^ = 0. 

Obviously, here lim^oo <p(t) = 0, lim^oo V((p(t)) — 0, so that the scalar field also goes to some limit (fioo, either 



finite or infinite, and V(ipoc) = 0. If the limit is finite, \poo\ < +oo, then, just by inequality (3.4), dV (ip) / d(p\ v=Voa = 0, 
d 2 V(ip)/dtp 2 \ v=lPoa > 0; the system evolves towards one of its finite stable critical (fixed) points with the zero minimum 
of the potential. It is not difficult to show that the expansion in this case is not accelerating unless p'(0) < —1/3, 
which means p + 3p < at least for small enough densities; the scalar field thus cannot "outweigh" the "normal" 
matter under such conditions. 

However, if tp^ = ±oo, which requires \im v —>± 00 V(<p) = 0, the situation with the acceleration can be different 
depending on the details of behavior of the potential at infinity (the system in any case goes again to a fixed point, only 
corresponding to an infinite value of the field) . We describe here two classes of potentials allowing the acceleration. 

a) Potentials vanishing exponentially at infinity. For some a > 1, let V(ip) = a(3a — 1) exp (—2p/^/a) [1 + o(l)] as 
ip — > +oo, and this asymptotic formula may be differentiated. If at l arge times <p goes to infinity, then for t — > oo, 



p ~ v^ln(t), p + 3p - (1 + 3v)/t 3a( - 1+ ^ with v = p'(0), and, by (gj), R/R ~ a{a - l)/t 2 if 3a(l + v) > 2, which is 
true for 1 + 3v > 0, i. e., exactly when the strong energy condition p + 3p > holds for small densities. The expansion 
in this case accelerates. 
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b) Potentials vanishing as a power at infinity. For some 6 > 0, let V{(p) = [4/(6 + i)] 2 (\/b / ip) b [1 + o(l)] as 
ip — > +oo, and this asymptotic formula may be differentiated. If at large times ^ goe s to infinity , then for t — > oo, 
<p ~ Vbt 2 /( b+4 \ p + 3p ~ (1 + 3i/) exp [-3(1 + i/)i 4 /< 6 + 4 )] with i/ = p'(0), and, by @, R/R ~ [4/(6 + 4)] 2 £- 2 /( b + 4 ), 
so that the expansion accelerates. 

These two types of potentials were the ones studied in |Ij - p] . 

To make the picture complete, we consider now the (not that interesting) exponential regime of expansion, i. e., 
solutions with p^ > 0. 

The large time behavior of the scalar field is described in this case by the damped anharmonic oscillator equation 
( p. 15 ). According to the Poincare-Bendixson theory of autonomous dynamical systems on the plane (e. g. ||, 2.8), 



every bounded non-closed phase trajectory has a limit trajectory which is cither a critical point or an isolated closed 



orbit (limit cycle, periodic solution). However, equation (0.15) has no periodic solutions at all. Indeed, it can be 
rewritten as 

p v = -2/30 2 , 

and the assumption that ip(t) is periodic with some period T > implies 

p v dt = Q= -2/5 / ip 2 dt, 
Jo 

that is, ip(t) = 0, because j3 > 0. The derivative of the scalar field is always bounded, <p 2 (t) < pt(/o) for t > to. 
Therefore if the scalar field is bounded, then t he so lution goes to a critical point, lirnt_»oo <p(t) = 0, limt_>oo <p(t) = 
<Poo, |<A»| < +oo, Viifoe) = p^? > 0, and, by (|0.15|) , dV(ip)/dip\ v=Voo = 0, d 2 V((p)/dip 2 \ v=Vaa > 0. Again the field, 



and with it the whole system, evolves towards a finite stable critical point which now corresponds to a non-zero 
minimum of the potential. This is the only possibility for a system with a growing potential, Iuiu-^oq V((p) = +oo, 
because under this condition all trajectories are obviously bounded. The system in such a case is dissipative, and its 
limit behavior is conveniently described in terms of the connected global attractor consisting of all the critical points 
and (heteroclinic) trajectories connecting the unstable critical points with th e sta ble ones ( pi, H])- 



It remains thus to find out about the unbounded phase trajectories of ( 0.15 ), if they exist at all, that is, un- 
der the condition that V((p) is bounded at least at one of the infinities. As explained, the unbounded trajectory 
means that lim^oo ip(t) — ±oo; to be precise, let us speak about the positive infinity. If the limit of the potential, 
lim^-f+oo V(cp) — Voo < +oo exists, then clearly lim^oo ip(t) = 0, = Pj? , so that the whole system once again 
tends to a stable fixed point, with just an infinite limit value of the scalar field. 

The exceptions of evolution towards a critical point thus could possibly be provided only by the potentials which 
are bounded but have no limit at infinity (it looks rather peculiar from the physical point of view, unless they are 
periodic). Our survey of the results pertinent to this case so far has shown that one cannot completely rule out the 
solutions which undergo rather strange damped oscillations characterized by the following: lim^oo ip(t) = ±oo, (p(t) 
is bounded and has no limit (in particular, lim^oo ip(t) ^ 0), lim^oo p v (t) — p§? > 0, and 

(p 2 dt < +oo, t > to 



If these exist at all, the relevance of such regimes to the Universe appears to be questionable at the very least. 

Let us now give a list of the most significant facts regarding the large time behavior of the cosmological expansion 
we have so far established. 

1. The matter density decreases to zero, the total energy density decreases to a (nonnegative) constant, the scalar 
field energy density is non- increasing and goes to the same constant as the total density, always. 

2. Except possibly for such potentials V(ip) that are bounded at infinity but have no limit there, the dynamical 
system always evolves towards one of its critical (fixed) points with either finite or infinite limit value (poc of the scalar 
field; the time derivative of the field goes to zero in both cases; if the limit value of the field is finite, the potential 
has a minimum at it. 

3. The expansion regime at large times is either exponential or subexponential with R/R — > 0; in the (typical) case 
of evolution to a fixed point the former occurs if V^oo) > 0, the latter if V(<poo) = 0. 

4. The subexponential expansion may be accelerating, that is, the scalar field can dominate "normal" matter 
(p + 3p > 0), provided that the scalar field tends to infinity and the potential has a zero limit V((foc) = V(oo) = 
and certain asymptotic behavior there. 

Properties 1 - 3 of the cosmological expansions are hardly surprising, however, all these results are now firmly and 
unambiguously established under general conditions of a clear physical origin. 
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Given the potential, the choice between the described possibilities is made by the initial conditions, so for the same 
potential the final stage of the expansion might be different depending on what happened in the early Universe. Let 
us illustrate this by some examples. 

If V(ip) is a potential well, i. e., it has a single minimum at some ip — ip* and goes (however slowly!) to infinity 
on both sides of it, all the solutions, disregarding the initial conditions they stem from, tend to the single (stable) 
equilibrium at tp*; depending on whether V(ip*) — or not, it is either a subexponential regime or an exponential 
one. 

A much more sophisticated picture arises with, say, the potential plotted in the Figure. It has a zero minimum 
at ip = if2, a non-zero minimum at ip = ip±, two unstable equilibria (maxima) at tp = P1.3, tpi < <f2 < </?3 < </?4, 
goes to infinity as ip — > +00 but tends to zero when ip — > —00. Both exponential (ipoo = (p^) and subexponential 
(Voo = ¥2, —00) regimes are possible, depending on how the expansion starts; the latter regime with ip^ = —00 
can be accompanied by acceleration even if p + 3p > 0, provided that V(tp) has the right asymptotic behavior at 
ip — ► —00. The exponential growth occurs, for instance, when tpo > ip% and p v (io) < ^(^3)) subexponential expansion 
without acceleration happens for the initial conditions with ipi < ipo < p% and p v (to) < V(ipi) < V(ips) (ip — ► ^2); 
subexponential regime with possible acceleration takes place when tpo < ipi and p v (to) < V(tpi), when p — > —00. 

Our final remark is about the closed Universe, k = 1. The only thing we need to validate our results in this case is 
that the expansion never stops (and maybe turns into contraction), i. e., that R > for any t > to- By (0.6), we are 
speaking about solutions for which pt — 3/i? 2 is positive for all t > to as soon as it is positive at the initial moment 
t = t ; for any such "infinite expansion solution", our results are true. There are many conditions under which these 
solutions exist. For instance, if V(ip) > V m in > 0, then, since pr > V(<p) by definition, any solution starting with 
Ro > Rrnin = \/3/V m i n is an infinite expansion solution. The other condition which does not depend on the potential 
at all is derived from the fact that, by the same definition, px > 8irp. The energy conservation equation (0.10) shows 
that at large times p ~ R~ 3i - 1+ "\ v = p'(0). If p'(0) < —1/3, then —3(1 + v) < 2, and any solution with large enough 
"starting radius" Ro satisfies pt > 3/i? 2 and hence is an infinite expansion solution, etc.. 

Therefore it seems that acceleration of the Universe is a fairly general behavior. What is surprising in light of recent 
supernovae results B is that the period of acceleration is just starting around the present time. One would have to 
fine-tune the parameters of our potentials to account for this. This rather unfortunate situation would be resolved if 
we had an appropriate potential which came from more fundamental physics. 

A paper with similar results for the technically much more complicated case of scalar-tensor theories of gravity with 
non-minimal coupling of the scalar field, is now in preparation. 
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